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FISSION-GAS RELEASE FROM URANIUM NITRIDE 
AT HIGH FISSION-RATE DENSITY 

by Michael B. Weinstein, Thomas A. Kirchgessner, and Thomas N. Tambling 

Lewis Research Center 

SUMMARY 

A sweep gas facility has been used to measure the release rates of radioactive fis- 
sion gases from small UN specimens irradiated to 8-atom-percent fissions (burnup) at 
high fission-rate densities. 

Fuel temperatures ranged from 600 to 1780 K. The measured gas release rates 
have been correlated with the predictions of a preliminary fission-gas transport- release 
model which includes atomic diffusion, fission- enhanced diffusion, and direct recoil 
release. 

Radioactive- gas release rates were found to increase linearly with fission burnups 
between 1. 5 and 8 percent. 

Postirradiation photomicrographs show a high- temperature (>1500K), high- burnup 

(6 percent) specimen to have large amounts of interconnected grain boundary porosity. 

137 

Measured cesium- 137 (Cs ) releases from the high-temperature, high-burnup speci- 
mens can be approximately predicted by assuming that the gas is released to the grain 
boundaries. 

Postirradiation photomicrographs after 6-percent burnup indicate that some bubble 
migration took place during operation at 1550 K. 


INTRODUCTION 

Efficient design of fuel elements for space power reactors depends in part on under- 
standing the fission-gas release and transport processes in these power- producing fuel 
elements. As part of an effort to increase our knowledge of these processes, a program 
to measure fission- gas release rates from operating nuclear fuels was initiated at the 
NASA Lewis Research Center. This study has as a goal the development of a method to 
predict fission-gas release rates from uranium nitride (UN). The technique being used 



is to measure radioactive fission-gas release rates from operating UN fuel specimens 
with a sweep gas facility and then to correlate these measured release rates with the 
test variables (temperature, fission-rate density, etc. ) through a theoretical model. 

The preliminary fission- gas transport- release model used in this study considers 
gas transport through the fuel by atomic and fission- enhanced diffusion (ref. 1), cul- 
minating in release of gas atoms from the surface. In this model, gas is also released 
by the recoil of gaseous fission fragments through the surface. 

Carroll, in a pair of review papers (refs. 2 and 3), indicates that the fission- gas 
transport and release mechanisms are complex and, as yet, not completely understood. 
Possible fission- gas transport mechanisms include atomic diffusion, fission- enhanced 
diffusion, and gas bubble diffusion, all of which can be complicated by trapping and re- 
lease processes. 

Atoms diffusing through the fuel can be trapped by defects or by stationary gas 
bubbles. They can be released from these traps by defect annealing or by kinetic reso- 
lution from the bubbles. Similarly, bubbles diffusing under the influence of a driving 
force, such as a temperature gradient, may be trapped or delayed by solid fission prod- 
uct precipitates, dislocations, or grain boundaries. Release from these traps occurs 
when the driving force exceeds the retarding force. 

These transport- trapping processes are themselves the product of multiple causes. 
As an example, Carroll (ref. 2) points out that possible bubble motion mechanisms are 
surface diffusion, volume (or bulk) diffusion, and evaporation- condensation. Domination 
of one bubble diffusion mechanism over the others depends on the material properties, 
the fuel temperature, and the bubble size (Nichols, ref. 4). 

Fission-gas release mechanisms include direct recoil (the ejection of fission- 
produced atoms from the fuel), knockout (ref. 5) (the ejection of gas and fuel atoms 
from the fuel surface by recoiling fission fragments), as well as the release when dif- 
fusing gas atoms or bubbles reach a free surface. Free surfaces include internal 
cracks or grain boundaries open to the exterior surface of the fuel. 

The tests reported herein were conducted to study the effects of three major vari- 
ables on fission-gas release. The variables are temperature, fission-rate density, and 
burnup. 

Small-size, high- enrichment UN samples were irradiated at high fission-rate den- 
sities to achieve high burnups (over .6 percent) in a reasonable length of time. Fuel 
temperatures ranged from 600 to 1780 K. 

Based on the diffusion transport model chosen (atomic plus fission- enhanced diffu- 
sion with desorption and recoil release), equations were developed which relate meas- 
ured fission- gas release rates to the test variables. Unknown constants in these equa- 
tions were evaluated from the data. 
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TEST APPARATUS 


The irradiation tests described in this report were conducted using a sweep gas ; 
facility located at the NASA Plum Brook Reactor Facility. A concise description, of the 
sweep gas facility, as well as the design of the fuel pin and the irradiation capsule, is 
presented here. More detailed descriptions can be found in reference 6. 


Sweep Gas Facility 

The sweep gas facility is basically a helium flow system which provides a means of 
monitoring the fission-gas release from power-producing fuel samples. A simplified 
flow schematic of the sweep gas facility is shown in figure 1. Bottle- supplied helium is 



Figure 1. - Simplified schematic of sweep gas facility. 
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Figure 2. - Instrument lead system. 



piped through a gas dehumidifier-purifier, a flowmeter, and an irradiation capsule con- 
taining a fuel specimen. Gaseous fission products released from the fuel sample are 
swept out of the capsule by the helium through a copper wool filter (to remove fission- 
produced iodine) and past a sodium iodide (Nal) scintillation detector, where the activity 
is measured. Samples of the sweep gas can.be removed for more detailed analysis. 
After flowing past the detector, the gases are stored in a large tank for intermediate 
holding and radioactive decay before release to the reactor off-gas cleanup system. 

The portion of the sweep gas facility located inside the reactor tank consists of an 
instrument lead system (fig. 2) containing the irradiation capsule, the iodine filter, and 
the solenoid valve assembly, plus the capsule insertion device. The entire instrument 
lead system is replaced as a unit with each new fuel specimen. 


Irradiation Capsule 

The small- volume irradiation capsule presently being used is shown in figure 3. It 
is designed to irradiate small- diameter fuel cylinders contained in a vented pin. 



Figure 3. - Irradiation capsule. 


The fuel sample is contained in a cylindrical molybdenum - 5-percent-titanium - 
0. 08- per cent- zirconium (TZM) pin (fig. 4) which is vented at the upper end to permit the 
escape of fission gases into the helium sweep gas. 

The fuel pin is axially supported and centered in the capsule by means of thermo- 
couple sheathing tubes (molybdenum) attached at each end of the pin. A third thermo- 
couple is set in a slot machined in the TZM clad. This latter thermocouple senses the 
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rEnd thermocouple End thermocouple -j 



clad temperature midway along the pin. 

The pin- capsule heat- transfer characteristics depend on the size of the gas gap be- 
tween the pin and the copper conduction cylinder. Fuel temperature is primarily a 
function of fission power for a particular gap size. 


Fuel Description 

Each fuel pin irradiated in this program contained four UN pellets having the follow- 
ing nominal dimensions: 0. 635 centimeter (0. 25 in. ) in length by 0. 127 centimeter 
(0. 05 in. ) in outside diameter by 0. 0763 centimeter (0. 03 in. ) in inside diameter. The 
total fuel length in each pin was approximately 2. 54 centimeters (1. 0 in. ). The fuel 
pellets were fabricated at Oak Ridge National Laboratory by using the general procedure 
described in reference 7. Some modifications to these procedures were necessary be- 
cause of the extremely small pellet diameters desired. 

The steps in the fabrication procedure are as follows : 

poc 

(1) Fully enriched (~ 93-percent U oa ) UN powder is isostatically pressed into small 
rods. 

(2) The rods are sintered at 2570 K for 2 hours in 1-atmosphere nitrogen (Ng) to a 
density of 13. 3 g/cm 3 (93 percent of theoretical density). 

(3) The sintered rods are centerless ground to a diameter of 0. 127 centimeter 
(0. 05 in. ). 

(4) The 0. 0763 -centimeter (0. 03-in. ) center hole is made by electrical discharge 
machining. 

(5) After inspection, the pellets are vacuum annealed. 

Pellet size is shown relative to a common pin in figure 5, and the typical fuel 
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Figure 6. - Section of fuel pin before irradiation. X500. 


microstructure is presented in figure 6. The typical physical, chemical, and isotopic 
characteristics of a pellet are as follows: 


Fuel UN 

Nominal length, cm (in.) 0.635 (0.25) 

Nominal outside diameter, cm (in.) 0. 127 (0.05) 

Nominal inside diameter, cm (in. ) 0.0763 (0.03) 

Percent density, percent of theoretical (theoretical density, 

14.3 g/cm 3 ) 93 

Oxygen concentration, ppm 450 

os* 

Enrichment, percent of U 93 

U 236 content, percent 0. 001 

_Q 

Grain size, cm ~ 1.5x10 
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TEST PROCEDURES 


All the irradiation tests discussed in this report were conducted in the following 
manner: During each 14- to 18-day Plum Brook reactor cycle, the fuel sample was 
irradiated at from one to six different temperatures with sweep gas pressure and flow 
rate held constant at all times. 

The sweep gas flow rate was usually held at 0. 2 cm /sec, so that the travel time 
from the capsule to the scintillation detector was 1130 seconds. This 1130-second 
travel time allowed the very short-lived fission-gas isotopes time to decay to an activity 
level too low to influence the scintillation detector data for the krypton and xenon iso- 
topes studied (Kr 87 , Kr 88 , Kr 85m , and Xe 135 ). 

Each constant- temperature (±10 K) run lasted from 24 to 48 hours, which was 
enough time for the release rates of the krypton isotopes to reach a quasi- steady state 
(release rate apparently unchanging in ~24 hours). Scintillation detector data were 
usually obtained every 6 hours (or less), and samples of the sweep gas were obtained 
every 3 to 7 days. 

There was continuous (strip chart) recording of two of the fuel- pin- mounted thermo- 
couples, as well as of the sweep gas flow rate. Hourly scans of all data channels were 
carried out automatically by the Plum Brook EDLAS (Experiment Data Logging and 
Alarm System). 

At the end of the total irradiation period the in-tank capsule assembly was trans- 
ported to the Plum Brook Hot Laboratory Facility for postirradiation examination. In- 
cluded in this examination were 

(1) Counting of capsule -mounted aluminum- cobalt wires for thermal fluence deter- 

minations 

(2) Counting of capsule-mounted stainless- steel and nickel wires for fast fluence 

determinations 

(3) Leaching of all internal capsule surfaces to detect vented uranium, cesium- 137 

(Cs* 37 ), cerium- 141 (Ce*^), and cerium- 144 (Ce*^) 

(4) Neutron radiography of the fuel pin to detect possible gross fuel changes 

(5) Fuel metallography to determine fuel density decreases and structural changes 

ooc ooc 

(6) Determination of IT °/U ratios for burnup determination 

137 144 141 

(7) Fission product (Cs , Ce , Ce ) counting to determine gross fission prod- 

uct retention in the fuel 

GENERAL TEST HISTORIES 
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A concise description of the irradiation tests is presented in table I. Each of the 


TABLE I. - OUTLINE OF IRRADIATION TESTS 


Capsule 

Total irra- 
diation 
time, 
hr 

Total 

burnup, 

percent 

Primary test 
temperature, 
K 

Temperature range 
(side thermocouple), 
K 

Fission-rate density 
range, 

3 

fissions/cm -sec 

Number of 
sweep-gas 
samples 

121 

1294 

8.3 

1780 

1215 to 1789 

3.34X10 14 to 8.29X10 14 

9 

122 

1923 

6. 6 

1223 

593 to 1225 

0.71X10 14 to 3.64X10 14 

19 

123 

1894 

6.0 

1505 

579 to 1512 

0.49X10 14 to 4. 45xl0 14 

19 


three capsules was tested primarily at one fuel temperature, with shorter periods of 
time at other temperatures. Complete temperature- time plots are given in figures 7 to 
9. The vertical lines on these plots represent temperature changes (planned or un- 
planned). 

Thermal and fast fluences for each test are given in appendix A. 



Figure 7. - Side thermocouple temperature as function of irradiation time - capsule 121. 
(Small variations not included. ) 
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Indicated temperature, K Indicated temperature. 



RESULTS AND DISCUSSION 


Results 

Our test and data-taking procedures are discussed in the section TEST PROCE- 
DURES. Data reduction procedures and an extensive sampling of fission-gas release- 
rate data are presented in appendix B. 

The release-rate data obtained in this series of irradiation tests can be correlated 
with a semiempirical model of fission-gas release which is based on concepts found in 
the literature. Test results and the correlation with the model are presented in this 
section. 

Many of the results are primarily based on data obtained in the tests of capsule 123, 
although any conclusions are supported by data from the other tests (capsules 121 and 
122; see table I). 

Change in release rate with burnup . - The basic krypton- 88^ release rate data for- 
capsule 123 are presented in figure 10 (similar data for capsules 121 and 122 can be 
found in appendix B). From the data in figure 10 it appears that, after about 1- to 1. 5- 
percent burnup (3.2x10 to 4. 8x10 fissions/cm ), the release rate increases with 
burnup and that the rate of increase may be independent of test history. 

A relation between release rate and burnup can be developed from the data in the 
following iterative manner: First, assume an approximate linear relation between re- 
lease rate (or activity in the sweep gas) and burnup above 1. 5 percent. This assumed 
relation is then used to group the data points in figure 10 in groups of constant burnup. 
Then the groups of data points at constant burnup but different temperatures are analyz- 
ed to obtain an average relation between release rate and temperature (table II), and 
this relation is then used to correct all the data points to a constant temperature (975 K). 

The temperature correction is then applied to the original data points, all at differ- 
ent burnups; and the data are analyzed to determine the best linear fit in a least- squares 
sense between activity and burnup. If the least- squares relation is found to differ from 
the original guess, it is used to repeat the cycle. 

When analyzed in this manner, the data in figure 10 yielded the linear relation 

r =r Q (l + 0.49 b) (1) 

where the symbols are defined in appendix E. 

The best- fit linear relation is shown in figure 11. All the data above a burnup of 


^Krypton- 88 is used for most of our analyses because its 2. 4-MeV gamma ray is 
the easiest to monitor and count. 
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Activity of Kr ° 8 in sweep gas, disintegrations/cm *-sec 







/- +20 Percent 



.5 1.0 L 5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 

Burnup, percent 


Figure 11. - Linear relation between release rate and burnup: activity of krypton-88 corrected to 975 K as function of 
burnup. Best fit: r/r 0 * (1 + 0.49 b). 
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Figure 12. - On-line activity of krypton-88 as function of burnup - capsule 122. Nominal tempera- 
ture, 1220 K. 



1. 5 percent in figure 10 are included. Almost all the data fall between limits of ±20 per- 
cent, and no consistent deviations with temperature were observed. 

Data from capsule 122 at 1220 K also agree with equation (1) as shown in figure 12. 

The low-burnup data (b < 2 percent) for all three capsules (figs. 10, 12, and appen- 
dix B) indicate a drop in release rate with increasing burnup to some minimum value 
prior to the beginning of the linear increase just discussed. A decreasing release rate 
has been observed by others (ref. 3, p. 566) and is generally explained as being caused 
by a decrease in surface area (refs. 3, 5, and 8) caused by surface smoothing or by 
sintering. We have also observed this effect during continuous irradiation at 800 K 

o 

(fig. 13) ; so it is not caused only by high- temperature operation. 



Irradiation time, hr 


Figure 13. - Decrease in krypton-88 release rate with time from start of irradiation. Temperature of capsule 1?4, 
~800 K. 


Though observed, this decrease in release rate with increasing burnup during low- 
burnup operation will not be included in our correlation because the complete dependence 
on temperature, fission-rate density, and burnup is unknown. 

Change in release rate with fission- rate density and temperature . - Based on the 
concepts of direct recoil of fission fragments, fission- enhanced diffusion (ref. 1), and 
atomic diffusion, an equation relating fission- gas release rates to fuel temperature and 
fission- rate density is 


dYSLf 

4T 


+ fYS 


' D fo f + D o e 


■Q/RT 


, 1/2 


( 2 ) 


This equation was derived as equation (C8) in appendix C. 


2 


Data analysis for capsule 124 is not complete. 
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The complete representation of the release- rate correlation combines equation (2) 
with the linear burnup relation (eq. (1)) to give equation (3). 


r = (1 + 0. 49 b) 


dYS M + fYs/M 


: + D. 


-q/rt\ 1//2 


b > 1 . 5 percent 


( 3 ) 


88 

By analysis of the Kr release-rate data for capsule 123 (fig. 10), values for the un- 
known constants in equation (3) were determined. The numerical form of equation (3) is 


r = (1 + 0.49 b) 


7.4X10 


r 4f y + fY /i.7xio" 29 f + 5.6x10“ V 156000 / 111 ^ 2 * 


(4) 


Values for d and Dj 0 were determined from the lower temperature (T s 1000 K) 
release-rate data. Next, the values for D Q and Q were determined from the higher 
temperature data; D Q and Q (representing atomic diffusion) are unimportant at the 
lower test temperatures. Equation (4) can be directly compared to the data presented in 
figure 10 by using the sweep- gas- activity - release-rate expression (eq. (Bl) of appen- 
dix B) 


r 


AQe XV//Q 


X 


(5) 


The solid lines drawn in figure 14 (which repeats the data in fig. 10) represent the values 
calculated from the semiempirical equation (4) and equation (5). The calculated lines 
shown in figure 14 agree closely with the measured data points. The mean deviation is 
within 11 percent for burnups over 1. 5 percent. 

When the constants determined from the capsule 123 data are used the calculated 
release rates are within ±14 percent of the measured values for capsule 121 (fig. 15) 
and are consistently higher than the measured values for capsule 122 (fig. 16). Data for 
these two capsules are given in appendix B. 

In appendix B it is pointed out that the gas sample activities have an estimated ac- 
curacy of ±10 percent, and that the calculated fission- rate densities are accurate to 
within an estimated ±15 to ±20 percent. Thus, the values found for the constants in 
equation (3) are no more accurate than these figures. 

We can compare the values found for some of the constants in equation (3) with val- 
ues reported in the literature. 
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Activity of Kr°° in sweep gas, dis integ rations/cm ^-sec 





Burnup, percent Burn up, percent 

Figure 15. - Comparison of empirical relation with data: activity of krypton-88 Figure 16. - Comparison of empirical relation with data.- activity of krypton-88 as 

as function of burnup - capsule 121. function of burnup - capsule 122. 





The value for the group of constants dSL/4 (direct recoil release) which best fits 

- 4 4 

our data is 7. 4x10 . If we use a recoil length in UN of 4. 8x10 centimeter and a 

geometric surface area of 0. 61 square centimeter, the fraction of recoils stopped by 

the helium at 330 K can be calculated (see appendix C) to be 

K = 7. 4x10”^ [ — — ) = 3. lx 10" 2 * 
\TSL/ 

Belle (ref. 9) reports that the fraction of recoils that are stopped by a gas gap between 
fuel and clad is given by 


l 

For our pin, if we use an h of approximately 0. 076 centimeter (i.e. , the specimen 

O 

inside diameter) and an l of approximately 2.3 centimeters, K ~ 2h/Z, or greater 
- 2 

than 6. 6x10 . The discrepancy between measured and calculated values for K is not 

actually a factor of 2. The true gas length traversed by a recoil fragment for our pin 
would be somewhat less than the 0. 076 centimeter used for h. 

Brinkman (as reported by Nelson in ref. 1) found that the fission- enhanced self- 

Q 1 

diffusion coefficient in uranium could be approximated by ~ 1. 5x10 f. Nelson 

states that this value could represent enhanced gas atom diffusion to within an order of 

-29 

magnitude. We have obtained a value of 1. 7x10 for the combination of terms 
2 

D f0 S , where S is the surface area of the fuel. Later, in our analysis of total gas re- 
lease from the fuel, we will show that - after assuming complete interconnection of 

on o 

grain boundary porosity - our value for D f would be approximately 1.3x10 cm /sec. 

QCm 

Isotopic variation in release rates . - On-line activity measurements for Kr° , 

87 185 8ft 

Kr , and Xe were made in addition to those reported for Kr . Analysis of sweep 

133 

gas samples added data for Xe . The release rates for all the measured isotopes 

Oft 

.showed the same increase with burnup as found for Kr , that is, v/r n = (1 + 0. 49 b), 

88 ° 

although the data analysis was not as detailed as for Kr . 

Another important observation is that the relative isotopic release rates do not ap- 
preciably shift with changes in temperature and fission- rate density. This is shown in 


2 Estimated for helium at 330 K and 6. 9xl0 5 N/m 2 (100 psia) from measurements in 

ref. 10 and density- temperature corrections in ref. 11. 
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Activity of Kr°° in sweep gas, disintegrations/cm 3 -sec 


Figure 17, - Isotopic activity ratios as function of krypton-88 activity. Temperature varies from 600 to 1550 K. De- 
cay constant, sec’ 1 : xenon-135, 2.09xl0’ 5 ; krypton-87, 1.48xl0 -4 . 

figure 17 where the activity ratio A ioc/A Qn is plotted as a function of krypton- 88 4 

Xe L ™ Kr° 

sweep gas activity, which increases with temperature, fission- rate density, and burnup. 

This observation seems to imply a single controlling fission- gas transport mechanism. 

The isotopic variations in release rates allow us to check equation (3), which pre- 

1/2 

diets an approximate variation of the form r a (1/x) ' . Combining this relation with 
equation (5) for isotopes 1 and 2 gives us the convenient form 

a 2 y 2\ x 2 / e - x 2 V ^ 

In table m we present the calculated and measured activity ratios for various iso- 
tope pairs. The results indicate that the variation with decay constant is roughly as 
predicted by equation (3). 

Postirradiation examination results . - Postirradiation fuel microstructures for 
capsules 123 and 122 are shown in figures 18 to 22. Comparison of figure 20 (capsule 
123, x500) and figure 22 (capsule 122, x500) with the preirradiation microstructure 
(fig. 6) shows that additional porosity has developed. In figure 22 this additional poros- 
ity is in the form of relatively few, small, intragranular voids, approximately 0. 5 to 
1 micrometer in diameter, uniformly dispersed throughout the fuel except near the in- 
side radial surface, where there is a high void density. Figure 20 (capsule 123) shows 
more and larger (0. 5 to 2 \im) intragranular voids preferentially located on the higher 
temperature side of each grain (side of grain facing right or toward inside surface). Fig- 

4 

Krypton- 88 activity was chosen for the abscissa in fig. 16 simply because the data 
analysis for this isotope was the most detailed. 
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TABLE HI. - COMPARISON OF CALCULATED AND 


MEASURED ISOTOPIC ACTIVITY RATIOS 


Isotope pair 

Calculated activity ratio, 
Aj/A 2 , if mechanism is - 

Measured a 

activity 

Tatio 

1 

2 

r oc (1/A) 1 / 2 

r oc (1/A) 

r * f (A) 

Kr 85m 

Kr 87 

0.487 

0.905 

0. 262 

0.584 

Kr 85m 

*-< 

oc 

oc 

.331 

.418 

.262 

.371 

00 

00 

U 

a 

t* 

OO 

Sh 

a 

1.48 

2. 17 

1. 01 

1.57 

b Xe !35 

Kr 87 

1.91 

5.07 

.717 

1.85 


a Average of 17 gas samples for capsule 123. 

k i nc 

°Neglects neutron absorption in Xe A 0 and iodine decay to 


xenon outside the fuel. 



Figure 18. - Section of fuel pin - capsule 123, X50. 





Figure 20. - Section of fuel pin - capsule 123, X500. 
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Figure 22. - Section of fuel pin - capsule 122, X500. 






ure 20 also shows interconnected intergranular porosity oriented perpendicularly to the 
direction of the temperature gradient. 

Figure 19 (capsule 123, x250) and figure 21 (capsule 122, x250) show the entire fuel 
thickness for each capsule and reinforce the previous observations about void number, 
size, and location. 

There has been some fuel swelling observed. The specimen swelling in capsule 123 

is a representative case and is discussed here. At assembly, the fuel- clad interface 

-3 -3 

gap varied from about 0. 6x10 to 1. 0x10 centimeter. Figure 18 shows the postirra- 
diation interface gap to be less than one- half the original size. The fuel inside diameter 
does not appear to have changed. The calculated volumetric increase (assuming iso- 
tropic fuel swelling) is about 3 to 6 percent. It is quite possible that the fuel was con- 
strained by the heavy walled clad. The small postirradiation interface gap seen in fig- 
ure 18 is what would be expected had the fuel swelled during irradiation, made contact 
with the constraining clad, and then separated on cooldown because of the unequal con- 
traction of the two materials. 

As mentioned in the section TEST PROCEDURES, the postirradiation examinations 

include measurements to determine vented uranium and gross fission product retention 

137 

and release (as indicated by Cs ). Results for the three irradiations are presented 
in table IV. 

The boiling point of cesium (943 K) indicates that cesium should have been gaseous 
during all our tests. Thus, it should be possible to calculate the long-lived cesium re- 
lease by extending our model to include stable gases. 


TABLE IV. - FUEL AND FISSION PRODUCT LOSS FROM FUEL SPECIMENS 



Capsule 121 

Capsule 122 

Capsule 123 

Primary test temperature, K 

1780 

1220 

1550 

Vented uranium 3, (in nitric acid 

0.3 

<0.2 


leach solution), mg 
137 

Average Cs retained, percent 

22 

(not detected) 
99 b 

55 c 

1 3 7 

Cs lost (in leach solution), percent 

65 

3.4 

40. 5 d 

141 

Average Ce retained, percent 

93 

96 


141 

Ce lost (in leach solution), percent 

0.2 

<1 

<0.1 


a Initial fuel weight for each capsule was ~0. 24 g. 

°A second sample showed 185-percent retention. 

C A third sample showed 15- percent retention. 

Efficiency of nitric acid leach not known; Cs 1 ^ losses are at least as high 
as given. 
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Equation (7) (which was derived as eq. (D26) in appendix D) gives the total stable 
gas release-birth ratio 


where 


_T _ KL + j _ _8_ 
B _ 4y 9 2 



1 - e 


-(2n+l)^a 


(2n + 1 ) a. 


n=0 


a A (D f + D)tt 


2 _t_ 
4 ?2 


( 7 ) 


Before this equation can be used to calculate r^,/B, we must obtain an estimate of the 
total surface area at the end of the irradiation period. The constants previously deter- 
mined in order to fit equation (3) to the radioactive- gas release rate data include com- 

2 2 

binations of terms like D^S and D Q S , where S and D Q are not separable. 

If we assume complete interconnection of the grain boundary porosity in the circum- 
ferential direction, the true surface area is approximately 14 square centimeters; and 
the correct length to use in equation (7) is the half-grain thickness (~7.5xl0~ 4 cm). 

Calculated Cs releases for capsules 121, 122, and 123 are compared to the 
measured releases in table V. The calculated releases do not include recoils (<1 per- 

TABLE V. - COMPARISON OF CALCULATED AND 

MEASURED CESIUM- 137 RELEASES 


Capsule 

Tempera- 

ture, 

K 

Time, 

t, 

hr 

Release 

parameter, 

a 

Ratio of total release to birth, r T /B 

Calculated 

Measured 

121 

1785 

460 

0. 859 




1780 

80 

.089 




1550 

170 

.066 




1225 

585 

.013 






= 1. 027 

0.48 

0.65 to 0.78 

122 

1225 

1375 

0.0315 

0.085 

0. 01 to 0. 03 

123 

1540 

680 

0.242 




1385 

180 

. 017 




1225 

215 

.004 






Y<* = 0.263 

0.25 

0.40 to 0.45 
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cent) nor any additional burst releases caused by abrupt temperature changes. 

137 

The discrepancy between the measured and calculated total Cs release for cap- 
sule 122 was expected since figure 22 does not show the large-scale interconnection of 
grain boundary porosity of the fuel in capsule 123 (fig. 20). 

While this analysis is very approximate in nature, it does indicate that the large 
fission- gas releases measured in these tests can be best explained by combining re- 
leases from individual grains with a gross interconnection of grain boundary porosity. 


Discussion 

Based on a simple recoil, atomic-diffusion, fission- enhanced- diffusion model, a 
fission- gas release- rate equation (eq. (3)) was derived and used to correlate the steady- 
state fission- gas release- rate data obtained with the sweep gas facility. Constants in 
the release-rate equation were determined from the test data taken with capsule 123 
(fig. 10 or^l4), so that for this capsule the mean deviation is within 11 percent (for 
burnups over 1. 5 percent). When the same constants are used, the calculated release 
rates are within ±14 percent of the measured values for capsule 121 and are consistently 
higher than the measured values for capsule 122, 

Both capsules 121 and 123 were irradiated to high fission burnups (8.3 percent and 
6. 0 percent, respectively) at temperatures over 1500 K. Capsule 122 was irradiated at 
a maximum temperature of 1225 K. As is shown by postirradiation photomicrographs 
(figs. 20 and 22), the UN irradiated in capsule 123 (fig. 20) contains large, 0. 5 to 2 ju m, 
gas bubbles plus interconnected grain boundary porosity, while the UN irradiated in 
capsule 122 (fig. 22) contains smaller bubbles with little or no grain boundary porosity. 

The gas transport- release model can also be extended to explain the high release of 
137 

Cs from the fuel if it is assumed that the fission gas is released first to the grain 
boundaries and then from the fuel. The presence of interconnected grain boundary 
porosity in the photomicrographs of the fuel in capsule 123 (figs. 18 to 20) supports this 
view. 

Recently reported gas release data for (Uq qPuq 2 )^ irradiations (ref. 12) has been 
interpreted as being related to surface- connected porosity. This interpretation has 
also been used to explain data obtained in UN irradiations (ref. 12). 

The (Uq glhiQ 2)^ data re P or ted in reference 12 show a steep increase in gas re- 
lease at about 1-percent burnup. At this point, pore interconnection is taking place. 
Test data for Kr^ release at 1550 K for capsule 123 (fig. 10) show a large increase at 
about 1-percent burnup - at which point, pore interconnection may be taking place. 

Our gas transport- re lease model does not include a description of directed bubble 
migration. Evidence for the existence of this mechanism was presented in figures 20 
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and 22. Consideration of this process as one which contributes to fission-gas release 
is consistent with our data. 

Analysis of high- temperature- gas release rates shows the exponential temperature 
dependence in eq. (2)), which we assumed to be caused by gas diffusion in the 

fuel. Nichols (ref. 4) shows that pore migration up a temperature gradient is dependent 
on temperature in approximately the same manner if the controlling mechanism is sur- 
face diffusion, volume diffusion, or evaporation- condensation. 

The UN fuel specimens irradiated in this program were operated with large (for UN) 
temperature gradients of as much as 500 K/cm. These gradients are a driving force 
for pore migration. 

A dynamic gas resolution process may also be important in describing fission- gas 
transport through the fuel. Gas permanently trapped in bubbles moving up a thermal 
gradient would be expected to attain a concentration varying approximately as 1/A. If 
no interchange existed between the gas held in the pores and the gas diffusing through 
the fuel, we could expect to observe a definite shift in relative isotopic release rates as 
the fuel temperature increased and the release controlling mechanism changed from 
atomic- plus- fission- enhanced diffusion to bubble diffusion. 

Our results, such as those shown in figure 17, indicate no shift in relative isotopic 
release rates when the fuel temperature changes from ^SOO K to 1550 K. 

This view of fission- gas transport in the fuel is consistent with ideas presented in 
many recent papers. For example, Carroll, in his review article (ref. 2), concludes 
that fission gas is in "constant turmoil" within fissioning fuel specimens. Gas resolu- 
tion and pore migration mechanisms can both be important in explaining fission- gas 
transport within fissioning fuel. 

We have noted that the high thermal gradients present in our fuel specimens would 
act as a driving force for pore migration. Similarly, our high fission- rate densities 
would increase the rate of gas resolution. 

The preliminary fission- gas transport- release model utilized in this report to pre- 
137 

diet the total Cs releases from high-burnup (>6 percent), high-temperature 

(>1500 K) UN specimens should not be directly applied to UN fuel pins operating at low- 

burnup and/or low- temperature conditions. A basic assumption in our utilization of 

137 

equation (7) to calculate Cs releases was that all the fission gas released from indi- 
vidual UN fuel grains was ultimately released from the specimen. Photomicrographs 
of specimen 123 operated to 6-percent burnup at 1505 K (figs. 19 and 20) show the 
presence of interconnected intergranular porosity. Gas released from individual grains 
can easily get to the exterior fuel surface. Without the formation of interconnected 
porosity, gas released to interior grain boundaries may be trapped and not released 
from the fuel pin. 

It is possible that low-burnup, low- temperature, fission- gas release can be de- 
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scribed by using our equations (such as eq. (7)) but considering only those grains which 
touch the geometric surface of the fuel. As the fuel burnup increases, more gas is 
released to the grain boundaries and finally some or all of this gas can escape through 
interconnected intergranular porosity. 


SUMMARY OF RESULTS 

This section contains a basic summary of experimental fission-gas release results 
obtained from small uranium nitride specimens irradiated to 8- percent burnup at high 
fission- rate densities: 

85m 87 88 

1. Release rates of the three krypton isotopes examined (Kr , Kr , and Kr ) 
decreased with increasing burnup during the low- burnup (0 to 1 percent) period of our 
irradiation tests. This decrease was observed at fuel temperatures from 800 to 1500 K. 

2. After a fuel burnup of approximately 1. 5 percent, fission- gas release rates in- 
creased linearly with burnup, independent of test history, to over 8-percent burnup 
(upper test limit). 

3. Postirradiation photographs of the capsule 123 fuel microstructure show grain 
boundary porosity oriented perpendicularly to the direction of the temperature gradient, 
as well as intragranular porosity preferentially located on the higher temperature side 
of each grain. 

4. The release rates of one isotope relative to another did not change appreciably 
with any test conditions. 

5. The variation of krypton release rate r with fission- rate density, temperature, 
and total burnup (b > 1.5 percent) for capsules 121 and 123 can be described by 

-Sg-iseooo/RTX 1 / 2 ! 


where f is fission-rate density, Y is fractional fission yield, T is temperature, R 
is the gas constant, and A is the decay constant. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1972, 

503-05. 


r = (1 + 0. 49 b) 


7.4xlO~ 4 fY + fy / 1.7xiQ~ 29 f + 5.6x10 
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APPENDIX A 


NUCLEAR ENVIRONMENT 


In order to more fully characterize the fuel specimen test conditions, the nuclear 
environment for each of the three specimens is outlined in table VI. The tabulated 


TABLE VI. - NUCLEAR ENVIRONMENT FOR SPECIMENS 


Capsule 

Maximum thermal 
neutron flux 
in fuel, 

neutrons/cm 2 - sec 

Maximum gamma 
heating rate 
in fuel, 

w/g 

Thermal neutron 
fluence, 
neutrons/cm 2 

Fast neutron fluence 
(>0. 1 MeV), 
neutrons/cm 2 

121 

5xl0 13 

3 

2xl0 20 

10 19 

122 

2 

2 

1 

10 19 

123 

3 

3 

1 

(Not measured) 


thermal- neutron- flux and gamma- heating- rate information is based on the specimen 
location in the reactor test hole during irradiation (with corrections applied depending 
on the reactor operating power and control rod bank height). Support tests were per- 
formed in the Plum Brook Mockup Reactor to determine the neutron-flux and gamma- 
heating-rate profiles along the length of the reactor test hole. A dummy capsule was 
used to obtain the information required to establish the capsule- specimen flux perturba- 
tion factors. 

Each test specimen was subjected to several thermal flux levels in order to achieve 
the fission heat generation rates required to maintain the several desired test tempera- 
tures. The tabulated values for thermal flux and gamma heating rate are the highest to 
which each specimen was subjected. 

The tabulated thermal and fast neutron fluence values were determined through the 
analysis of dosimeter wires located within and at the outside surface of each capsule. 
Aluminum- cobalt wires were used for determining the thermal fluence. Nickel and 
stainless- steel (300 series) wires were used for the fast fluence determinations. The 
thermal fluence dosimeters were located an appreciable distance from the fuel specimen. 
Hence, corrections were applied to the measured values in order to determine the ap- 
propriate thermal fluence values for the fuel specimens. The required correction fac- 
tors were determined in the Mockup Reactor tests. 
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APPENDIX B 


DATA AND DATA REDUCTION 

Of the variables of interest (release rate, temperature, fission-rate density, and 
burnup) only one, the temperature, was obtained from direct measurements. Fission- 
gas release rates were calculated from measured activities in the sweep gas, either 
from samples counted by the Plum Brook Radiochemistry Section or from on-line scin- 
tillation detector data. Absolute fission- rate densities were calculated from measure- 
ments of relative thermal flux as a function of position in test hole (relative fission- rate 
densities) combined with the total burnup obtained in the postirradiation examinations. 
Burnup could be calculated at any time during the irradiation by knowing total burnup 
and test history. 

An understanding of the accuracy to which the variables are known can be obtained 
by following through the sequence of steps in the data reduction procedure, as is done 
in this section. 


Sweep Gas Activities and Fission-Gas Release Rates 

Fission-gas activities in the helium sweep gas were determined by gamma-ray 
counting techniques, by using either an on-line scintillation detector or a Ge(Li) detec- 
tor for the gas samples. Analysis of the gas samples was estimated to give activities 
accurate to ±10 percent for the isotopes Kr 85m , Kr 87 , Kr 88 , and Xe 135 . An unknown 
error is introduced here; literature values for absolute gamma- ray intensities (gammas 

of a particular energy per disintegration) must be used. The values used are shown in 

87 

table VII, and particular note of the range of values reported for Kr should be made. 

TABLE VII. - GAMMA- RAY INTENTISY DATA 


Isotope 

Gamma- ray 
energy, 
MeV 

Gammas per 100 
disintegrations 

Reference 

Kr 85m 

0.150 

75.2 

13 

Kr 87 

.403 

a 59. 6 

14 



50 

15 



~68 

16 

CO 

CO 

t* 

2.392 

37.8 

17 

Xe 135 

.250 

97 

18 


a Value used in data analysis. 
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TABLE VIII. - FISSION-GAS ACTIVITY IN HELIUM SWEEP GAS - CAPSULE 121 GAS SAMPLES 3 


Fuel maximum 
temperature, 
K 

Fission- rate density, 
fissions/cm^-sec 

Fuel 

burnup, 

percent 

Isotope activity, 

3 

disintegrations per cm per second 

Kr 85m 

Kr 87 

5 

00 

oo 

Xe 135 

1236 

3. 43xl0 14 

0.54 

30.2X10 3 

55.9X10 3 

63.4X10 3 

69.8X10 3 

1240 

3.47 

.82 

29.7 

52.7 

59.7 

64.5 

1550 

5.77 

1.7 

262 

398 

490 

726 

1765 

7.97 

5.06 

1320 

1780 

2560 

3650 

1765 

7.97 

6.65 

2100 

2280 

3690 

5810 

1215 

3.30 

8.3 ' ; 

95.9 

184 

183 

227 


a Some samples obtained during each capsule test grossly disagreed with on-line data, probably 
because of leaks or incomplete filling. These have not been tabulated. 


TABLE IX. - FISSION-GAS ACTIVITY IN HELIUM SWEEP GAS - 
CAPSULE 122 GAS SAMPLES 


Fuel maximum 
temperature, 
K 

Fission- rate density, 

3 

fissions/cm -sec 

Fuel 

burnup, 

percent 

Isotope activity, 

3 

disintegrations per cm per second 

Kr 85m 

W 

00 

-3 

CO 

00 

. Xe 135 

1222 

3.64X10 14 

0.49 

24. lxlO 3 

38.8X10 3 

60. 5xl0 3 

50. lxlO 3 

1209 

3.43 

1.48 

20.5 

33.3 

45.5 

45.3 

1223 

3.64 

2.64 

30.6 

50.7 

73.5 

68.3 

1000 

2.33 

4.31 

15.8 

24.5 

38.5 

41.2 

1225 

3.64 

5.38 

40.3 

65.5 

104 

106 

800 

1.44 

5.81 

8.34 

15.9 ‘ 

25.7 

25.2 

1000 

2.33 

5.93 

22.4 

35.3 

57.2 

58 

598 

.71 

6.16 

3.20 

5.35 

8.00 

8 . 04 

805 

1.45 

6.27 

10. 1 

16.5 

26. 1 

26.3 

1.227 

3.64 

6.59 

44.3 

71.2 

118 

121 


The fission- gas activities in the sweep gas as determined from gas samples for 
each of the three specimens tested are presented in tables VIII to X along with other 
pertinent test data. A discussion of how we obtained the values for fission- rate density 
is presented later in this appendix. 

A typical scintillation detector spectrum is shown in figure 23. Integration under 
the peaks corresponding to the different isotopes yielded relative activities which could 
be converted to absolute activities by comparison with gas samples. Activities deter- 
mined by each method are compared in table XI. 

On-line activity data for each of the capsules are plotted as a function of burnup in 
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Fuel maximum 
temperature, 

Fission- rate density, 
fissions/cm 3 -sec 

Fuel 

burnup, 

Isotope activity, 

3 

disintegrations per cm per second 

K 


percent 

Kr 85m 

Kr 87 

* 

co 

00 

Xe 135 

1540 

4. 47xl0 14 

0.75 

41. 1X10 3 

71X10 3 

104X10 3 

131X10 3 

1536 

4.46 

1.42 

79.7 

144 

206 

237 

800 

1. 09 

1.54 

4.36 

6.79 

10.7 

13.3 

1534 

4. 45 

2.59 

185 

326 

504 

556 

1543 

4. 49 

3.45 

256 

426 

685 

993 

580 

.48 

3.52 

1.84 

3.85 

5.50 

5.97 

801 

1. 09 

3.63 

6.20 

10.8 

16.4 

19. 1 

992 

1.77 

3.72 

16.8 

26.6 

44. 1 

55.3 

591 

.50 

4.33 

2.49 

4. 85 

6.68 

8. 14 

812 

1. 12 

4.39 

8.24 

13.2 

22.0 

26.4 

1145 

2.43 

4.88 

39.9 

67.7 

109 

127 

995 

1.78 

5.04 

17.6 

29.8 

46.8 

51.5 

1214 

2.75 

5. 19 

36.9 

65.7 

108 

103 

1382 

3.61 

5.39 

122 

240 

351 

369 

1540 

4.47 

5.86 

357 

555 

939 

1200 





















TABLE XI. - Kr 88 ACTIVITY: ON-LINE COUNT 
RATE COMPARED TO ACTIVITIES 
MEASURED IN GAS SAMPLES 


[Data for capsule 123 !] 


Sample 3, 

Reported activity, 
disintegrations/cm -sec 

On-line count rate 
at 10 cm, 
counts/sec 

1 

1.04X10 5 

110 

2 

2.06 

250 

3 

1.07X10 4 

9.65 

4 

5.04X10 5 

570 

5 

6. 85 

700 

6 

5.50X10 3 

5.35 

7 

1.64X10 4 

14. 8 

8 

4.41 

35.0 

10 

6.68X10 3 ■ 

5.0 

11 

2.20X10 4 

17.0 

13 

1. 09x1 0 3 

89 

15 

4. 68X10 4 

42.0 

16 

1.08X10 5 

130 

17 

3.51 

350 

18 

9.39 

950 


Average ratio, On-line count rate = g.^xio' 4 
Sample activity 

Mean deviation from average ratio = ±12 percent 


a Unused samples show gross disagreement, probably 
because of leaks. 


figures 24 and 25. Nominal temperature is a parameter. Gas sample points have also 
been included. 

Sweep gas activities were converted to release rates by using the equation 


r 


AQe AV /Q 

A 


(Bl) 
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Absolute Fission-Rate Densities 


Absolute fission- rate densities corresponding to particular fuel-pin temperatures 
and capsule positions in the Plum Brook Reactor test hole RA-5 were calculated by 
combining estimates of relative thermal flux (assumed proportional to relative fission- 
rate density) with postirradiation measurements of total burnup. 

For particular capsule test positions we calculated the relative thermal flux by 
using the thermal flux profile obtained through Mockup Reactor measurements (fig. 26) 
and corrections for reactor shim- rod position and total power level. Plots of tempera- 
ture as a function of relative fission- rate density were then prepared for each reactor 
cycle by assuming that relative fission- rate density equaled relative thermal flux. 
Cycle-to- cycle variations in reactor characteristics caused by changes in core loading 
or experiments were eliminated by equating the relative fission- rate densities at one 
particular temperature. Figure 27 shows a plot of temperature as a function of fission- 
rate density for capsule 123. For each irradiation cycle, the relative fission-rate den- 
sity at 975 K was set to unity. 



Figure 27. - Temperature as function of relative fission-rate density with each cycle separate - capsule 123. 


3 This Mockup Reactor closely simulates the characteristics of the Plum Brook 
Reactor (PBR). 
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Absolute fission- rate densities were obtained by integrating the relative fission- rate 
densities over each pin’s irradiation history and then utilizing the measured total burnup 

O 

(fissions/cm ) to get the necessary conversion factors. Figure 28 shows the relation 
between temperature and fission-rate density for all three capsules. 

We have made no corrections for variations in fuel content during irradiation. For 
all our tests, any corrections would be only a few percent, or much less than the uncer- 
tainties inherent in this calculational procedure (15 to 20 percent). 
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APPENDIX C 


STEADY-STATE FISSION-GAS RELEASE RATE 


The general time- dependent equation for diffusional transport of fission gas in a 
one- dimensional rectangular coordinate system is 


— = (D. + D) — - Ac + fY 


at 


3y 


- (Cl) 


where it has been assumed that the diffusion coefficients Dj and D are constants and 
that the fission- rate density f is independent of position in the fuel. 

The steady- state form of equation (Cl) is obtained by setting 

— = 0 (C2) 

at 

And the two boundary conditions required to get a complete solution to the steady- state 
equation are 

— = 0 at y = 0 (C3) 

, dy 

and 


c = 0 at y = y 1 


where y^ is the distance across the fuel (fuel thickness). 

Descriptively, dc/dy = 0 at y = 0 represents a barrier to fission-gas transport, 
such as the fuel clad interface; and c = 0 at y = represents the conditions at the 
fuel- gas interface. 

The steady- state solution to equation (Cl) with the boundary conditions (C3) is 



(C4) 


where F and E are defined as 
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F A fY 
D f + D 


E 4 — - — 

Dj + D 

Now we assume that at the fuel surface (y = y^) fission- gas release is given by a 
recoil component added to the diffusional release.. The release equation is 

r = - (D. + D)S — (C5) 

4 dy 

y Yx 

where we have assumed that one-fourth of the recoil fragments originating within a dis- 
tance L of the surface leave the fuel. Combining equations (C4) and (C5) with the sim- 
plification 


V E 


*1 


» e 


-V E 




gives 


r 


,SSM + (ys 


(C6) 


Dj can be put into the form (ref. 1) 


where Dj q is a constant. D is usually presented as 

D = D 0 e-^ RT 

With these relations, equation (C6) becomes 


r - KYSLf 
4 


+ fYS 



(C7) 
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The constant K in equation (C7) is approximately proportional to the density of the 
gas into which the recoils are ejected (ref. 9, p. 652). Therefore, to describe the 
change in gas density with temperature, we can replace K by the term d/T, where d 
is a constant, to give the final form of the steady- state release- rate equation 


r ,dYSLf +fYS 

4T 


rD fo f + 


D e 'Q/ RT V /2 


(C8) 
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APPENDIX D 


TIME-DEPENDENT STABLE GAS RELEASE 


The starting point for this derivation is the one -dimensional transport equation (Cl) 


— = (D. + D) — - Ac + fY 

ay 2 


at 


(Cl) 


which for a stable gas (X = 0) becomes 


, 3 2 

— = (D f + D) — + fY 


at 


(Dl) 


sy 


Carslaw and Jaeger (ref. 19, pp. 29-30) state that the solution to an equation of this 
type with initial and boundary conditions 


c(y, 0) = M(y) 
c (surface) = N(y2) 




j 


(D2) 


can be expressed as 


c(y,t) = u(y) + w(y,t) 


(D3) 


where the functions u and w are themselves solutions for 


a 2 u -fY 
ay 2 D f + D 


(D4) 


and 


a 2 w l aw 

3y 2 D f + D 9t 


= 0 


(D5) 


with initial and boundary conditions 
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u (surface) = N(yr) 


w(y,0) = M(y) - u > 


w(surface) = 0 


(D6) 


We will first derive the solution for c(y, t) with an arbitrary initial concentration 
M(y) and boundary conditions 


c(surface) = c(y 2 ) = 0 = N(y 2 ) 


— = 0 at y = 0 
9y j 


(D7) 


Then this solution will be specialized to the cases 


M(y) = 0 


(D8) 


(zero initial concentration) and 


M(y) = c(y, r) (D9) 

where c(y, t) is the concentration found by using condition (D8) for a specific set of con- 
stants Dp D, and f at time r. Equations for stable gas release rates can be obtained 
for the special cases. 


Arbitrary Initial Concentration 


As given by equation (D3), the concentration c is 


c(y,t) = u(y) + w(y, t) 


(D3) 


where u(y), as defined by equation (D4) with boundary condition (D7), has the solution 


u(y) =■ 


fY 


2(D f + D) 


(A - y 2 ) 


(DIO) 
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Then w(y,t) is the solution to equation (D5) 


d 2 w _ 1 9w _ q 

Sy 2 D f + D 9t 


(D5) 


with initial and boundary conditions given by equations (D6) and (D7), which for this case 
are 


w(y, 0) = M(y) - 


2(D f + D) 
w(y 2 ) = 0 


(yf - y 2 )'! 




9w 

9y 


= 0 at y = 0 


(Dll) 


Equation (D5) may be solved by separation of variables. Let w(y, t) be given by 


w(y, t) = G(y)H(t) 


(D12) 


By differentiation, equation (D5) becomes 


H d 2 G G dH = 0 ^ 
dy 2 D f + D dt 


or 


d 2 G/dy 2 = dH/dt = _^2 
G H(D f + D) 


■ (D13) 


The general solutions to these equations are 


H =aje 


- (D f +D)j3 2 t ^ 


and 


G = a 2 sin j3y + a^ cos /3y 


(D14) 
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but since by condition (Dll) 



at y = 0 


a 2 = 0 and w(y, t) is then equal to 


- (Dj+D)/3 2 t 

w(y, t) = a 4 e A cos j3y 


where 


a 4 “ a l a 3 

Another boundary condition (Dll) is 

w(y 2 ) = 0 


or 


a 4 cos /3y 2 = 0 

The values of 0 for which this is true are given by 


0 


n 


2n + 1 n 

2 y2 


n = 0, 1,2. . . 


Thus, the solution (D15) becomes 


w(y,t) 


E -<D f+ D)$ 

V 

n=0 


COS 0 n y 


where the terms form an orthogonal system. The constants a n 
initial condition (Dll), which gives 


w(y, 0) = M(y) - 


fY 


( A. ■ y2 ) 


- Z a n C0S ^ 
n=0 


(D15) 


(D16) 


(D17) 


(D18) 

can be found by usir^ 


Y (D19) 


2(D f + D) 
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from which we obtain 


.1 

n 


72 r fY(y2 - y 2 \ 

M(y) ^ ^ 


2(D f + D) 


cos /3 n y dy 


f 2 COS 2 /3 n y dy 


(D20) 


Integrating equation (D20) gives 


n y 2 


f 


M(y) cos 0 n y dy 2fY - ~ — 

y 2 ^( D f + °) 


(D21) 


The concentration c(y, t) can then be obtained by combining equations (D3), (DIO), (D18), 
and (D21). 


c(y,t) = 


_2 

y2 


f 


M(y) cos By dy 2f ^ -1 ^ — 

n o 


y 2 ^( D f + D) 


■ (D f +D)/3^t 


cos /3 n y 


n=0 


fY / 2 2 


2(D f + D) 


(y 2 - y 2 ) (D22) 


Special Case - Zero Initial Concentration 


The initial condition for this special case is given by equation (D8) 


M(y) = 0 


(D8) 


Applying this condition to equation (D22) gives 
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oO 


c(y,t) =• 


fY 


2(D f + D) 


2 2 4 

y 2 -y - — 
?2 


_ -(D f +D)/3“t 
(-l) n e f n cos /3 n y 


n 


n=0 


which, if rearranged slightly, yields 


c(y,t) = 


fYy 2 

2(D f + D) 


l-z!- 32 
?2 n Z 


- (Dr+D)£Tt 

(-l) n cos/3 n ye f n 
(2n + l) 3 


n=0 


(D23) 


Equation (D23) is equivalent to the heat- transfer solution presented by Carslaw and 
Jaeger (ref. 19, p. 130). The release rate from the fuel is related to the concentration 
gradient at the fuel surface by equation (C5) 


r _ KYSLf 
4 


- (D, + D)S ^ 
1 dy 


*2 


(C5) 


where 


dc 



is obtained by differentiation of (D23). The release rate is 


r 


KYSLf + fY ■ s 
4 £ 


1 




- (D f +D)/3 3 t 
e 

(2n + l) 2 


(D24) 
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The total amount of gas released from t = 0 to t = t is found by direct integration 
of equation (D24): 


\ 

e -(2n + l) 2 a _ 
(2n+ l) 4 -^ 


where we have defined 



(D25) 


a # (D f + D)tt 


2 t 


This result can be expressed in the form 

Total release 
Total birth 

by dividing by B = fYygSt, which then gives 


r T . KL 
B ~4y 2 



j _ e -(2n+l) 2 Q! 
(2n + 1 ) 4 0 ! 


(D26) 


Equations (D24), in terms of r(t)/r(t = °°), and (D26), neglecting recoil, are plotted in 
terms of a in figure 29'. 
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Special Case - Step Change in Conditions 


This solution corresponds to a step change in fuel operating conditions (D^ - Dp 
D — D’, f — f’) after operating for some time r starting at an initial concentration of 
zero. 

For this case, the concentration M(y) is given by equation (D23). The concentra- 
tion c(y, t), where t is now the operating time after the step change at r, is found by 
combining equation (D22) with (D23), where D, Dj, and f have been replaced by D’, 
D' v and f': 


c(y,t) 



2 _ 

?2 



M(y) cos /3 n y dy - 


2f'Y(- l) n 

y 2 4< D f + D '> 


e 


cos 0 n y 


f'Y 


2(D| + D') 


- (y| - y 2 ) (D22) 
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M(y) =■ 


fY Y2 


2(D { + D) 


1- yI-32 
2 3 

y 2 77 


(-1) 1 cos /3.y e 


■ (D f +D)/3? r 


( 2 i + ir 


(D23) 


and the summation index has been changed to i. Integration yields equation (D27) 


c(y,t) = 



2fY(- l) n 


2 - 


-(D f+ D)/3*r 

1 - e 1 n 


2f'Y(- l) n 


y 2 ( D f + d)^ 


y 2 4< D f + D’) 


-(d:+d')A 
e cos /3 n y 


f’Y(y| - y 2 ) 

+ ^ (D27) 


2(D| + D’) 


The release rate is again obtained from equation (C5) 


r= KYSLf'- (D , + D , )g dc 
4 1 dy 


?2 


(C5) 


which yields 


r= _KJ^Lf' + s + 

A & 



n=0 


-\ 

S~ 

- (D»+D) (3 2 t 
1 - e 1 n 


2fY 

(D| + D’)S< 


_ 

y 2 (D f + D >4 

J 




2f'Y 


y 2 /3 n( D f + D ’) 


/n e 


-(D^+D’)/3 2 t 


(D28) 


50 


appendix e 


SYMBOLS 

Q 

A activity, disintegrations/cm -sec 

Aj activity of isotope 1, disintegrations/sec 

Ag activity of isotope 2, disintegrations/sec 

a constant 

a n constant coefficients 

B total birth (atoms formed by fissions), atoms 

b burnup, percent uranium atoms fissioned 

Q 

c concentration of fission products, atoms/cm 

o 

D diffusion coefficient, cm /sec 

D' diffusion coefficient for second time step 

Dj fission-enhanced-diffusion coefficient, cm /sec 

Dj fission-enhanced-diffusion coefficient for second time step 

D q diffusion constant, cm /sec 

Dj 0 fission-enhanced-diffusion constant, cm /sec 

d combination of terms, KT 

E combination of terms, X/CD^+D) 

F combination of terms, fY/(Dj + D) 

O 

f fission-rate density, fissions/cm -sec 

f' fission-rate density for second time step 

G a function of y, defined by eqs. (D12) and (D13) 

H a function of t, defined by eqs. (D12) and (D13) 

h gas gap thickness, cm 

K fraction of recoils stopped by gas after ejection 
L recoil length in fuel, cm 

l recoil length in gas, cm 

O 

M initial concentration in fuel, atoms/cm 
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2 

N surface concentration, atoms/cm 

Q energy of activation for diffusion, J/mole 

Q volumetric flow rate, cm /sec 

R gas constant, 8. 314 J/mole-K 

r release rate, atoms/sec 

r Q release rate at time equals zero, atoms/sec 

r,p total release, atoms 

o 

S surface area, cm 

T temperature, K 

t operating time, sec 

u a concentration function, defined by eq. (D4) 

3 

V sweep gas system flow volume, cm 

w a concentration function, defined by eq. (D5) 

Y fractional fission yield, atoms/fission 

y thickness, cm 

y^ slab thickness, cm 

y 2 slab half thickness, cm 

a release parameter, (D f + D)?r t/4y 2 

/3 a constant 

0 eigenvalues, defined by eq. (D17) 

X decay constant, sec” 1 

r operating time, total time for first set of conditions, sec 
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